THE UNIVIFSTTY OF MIAMI 



The Dyn?^nic Sti'v.cl:ure of a Wind -Induced Eddy 



BY 

Thomas Stephen O'Keefe, Jr. 



A THESIS 

Submitted to the Faculty 
of the University of Miami 
in partial iulfiilment of the requirements for 
tne degree of Master of Science 



Cor'il Gables, Florida 
August 1971 




.0 363 



.AL POSTGRADUATE SCHOaU 
CiJalF. 93940^>‘ 



O’KEEFE, THOMAS STEPHEN, JR, 



(M.S., Physical Oceanography) 



The Dynamic Structure of a Wind-Induced Eddy , (August 1971), 

Abstract of a Master’s Thesis at the University of Miami. Thesis 
supervised by Professor Claes Rooth. 

The response of the ocean to an axisymmetric cyclone is 
investigated using a model based on the conservation of potential 
vorticity in the interior. The results, which confirm the conservation 
potential vorticity, are used to obtain limiting values of the 
vertical eddy diffusivity and an effective Ekman depth. Some aspects 
of the response as a function of the horizontal wavenumber are al^o 
investigated. It is found that high wavenumber motion does not penetrate 
as deeply as low wavenumber motion. 



THE UNIVERSITY OF MIAMI 



The Dynamic Structure of a Wind -Induced Eddy 



BY 

Thomas Stephen O'Keefe, Jr. 



A THESIS 



Submitted to the Faculty 
of the University of Miami 
in partial fulfillment of the requirements for 
the degree of Master of Science 



Coral Gables, Florida 
August 1971 






, 6 > 



THE UNIVERSITY OF MIAMI 



A thesis submitted in partial fulfillment of 
the requirements for the degree of 
Master of Science 



Subject 

The Dynamic Structure of a Wind -Induced Eddy 



Thomas Stephen O'Keefe, Jr. 



ACKNOWLEDGEMENTS 



During the time I studied at the School of Marine and Atmospheric 
Science, I was fortunate to develop lasting friendships with many 
people. I have also experienced a perceptible change in attitude 
toward science and have developed a sympathy for the problems of the 
scientist, and a respect for his motivation. The answer to the 
quest ion**why** remains man*s greatest challenge. The friendships 
made, and the change in perspective will be the most lasting and 
beneficial effects of the last two years. 

I am indebted to the members of my thesis committee. Dr. Kamal 
Yacoub and Dr. Walter Duing. I am particularly grateful for the 
infinite patience of Dr. Claes Rooth, the committee chairman, who 
guided my efforts. I am thankful for the provocative advice of 
Dr. C.N.K^ Mooers who gave me a particularly valuable commodity, 
his time. 

Thanks are due Mr. Manuel Bascuas, who provided programming 
assistance, and Miss Linda Monaco, who typed the several drafts of 
the manuscript. 

Finally, I thank my wife, Kathy, for understanding. 

The financial support for this study was provided by the Office 
of Naval Research Contract No. N00014-67-A-0201-0013 , Project No. 

NR 083-060. 



iii 



TABLE OF CONTENTS 



Page 

LIST OF FIGURES v 

LIST OF APPENDICES vi 

I. INTRODUCTION 1 

II. METHODOLOGY 5 

III. DISCUSSION OF RESULTS 12 

A. Discussion of the Solution and 

Consistency of Results 12 

B. The Mixing Process and Momentum Transfer 18 

C. The Effect of Length Scale on Response 22 

D. Geostrophic Velocities 25 

IV. SUMMARY AND CONCLUSIONS 29 

APPENDIX 31 

LIST OF REFERENCES 39 



iv 



LIST OF FIGURES 



FIGURE Page 

1. The path of Hurricane Kara and stations of 

R/V PILLSBURY. Wind speed in knots is indicated 2 

2. The relation between density (p) , undisturbed 

depth (Zq) and stretching (§) 7 

3. Observed mass structure of ocean after passage of 
Hurricane Kara. Isopleths are values of 

sigma -t (C7|-) 11 

4. Comparison between the observed mass structure 
and the computed response after Hurricane Kara. 

Jo(kr) = 0 at 80 km 13 

5. Log-linear plot of temperature (®C) vs depth (m) 19 

6. Density difference (gm/cm^) plotted vs square of 
depth (cm^) for potential energy calculation. 

Station #24 is compared to the center of the 

computed curves ..... 21 

1 . Isopleths of stretching, §(m) , as a function of 

depth and length scale. Values of § are normalized 

to 100 meters pumping at 25 meters matching depth 24 

6. Comparison of numerically calculated geos trophic 
velocities (cm/sec) with velocities found using 
the dynamic method (observed) 27 

9. 0 = constant surface separated into horizontal and 

vertical components to illustrate coordinate 

transformation technique 33 



V 



LIST OF APPENDICES 



APPENDIX Page 

A. The Natural Coordinate System 31 

B. The Potential Vorticity Conservation Model 34 



I, INTRODUCTION 



In October of 1968, Tropical Storm Kara forced R/V PILLSBURY 
south from an operating area in the vicinity of latitude 30°25*N, 
longitude 76°52*W. While in the stages of maturation, the tropical 
storm displayed the curious behavior of slowing and crossing its own 
path several times. Kara eventually reached hurricane force, increased 
speed and moved on a generally northeast course. In transit to the 
operating area after passage of the storm, the research vessel 
experienced set to the right and then left, indicating that a hurricane- 
induced eddy had been crossed. Bathythermograph stations were taken 
in order to locate the associated dome-like density structure, and STD 
stations across the presumed dome were taken after it had been located. 
The track of Kara and the track of R/V PILLSBURY are shown in Figure 1. 
The hurricane track was obtained from the advisories issued by the 
National Hurricane Center. 

Kara may be considered to have been quasi-s tationary during the 
period 14-15 October 1968. Therefore it affords an opportunity to 
study some aspects of hurricane-induced dynamics through the upper 
part of the water column. As Stevenson (1966) points out, hurricanes 
are unique "laboratories" for the study of air-sea interaction and 
ocean dynamics. Although the significant atmospheric effects of a 
hurricane are spread over horizontal scales on the order of several 
hundred kilometers, the most interesting oceanic effects take place 
within an area which can be crossed in one day. Unfortunately, 
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Figure 1. The path of Hurricane Kara and stations of R/V PILLSBURY. Wind speed in knots is indicated. 
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investigations of the response to hurricanes depend by necessity on 
**cruises of opportunity,** Given the availability of research vessels, 
studies of this sort are further complicated by the unpredictability 
of the occurrence, duration and track of a storm. 

A study of low intensity storms such as Kara is particularly 
relevant because the wind velocities are not too different from those 
found in the winter storms in the North Atlantic Ocean. In particular, 
the data provide a means of testing assumptions of conservation of 
potential vorticity in the ocean. An analysis of those areas where 
the conservation equations are not applicable is important because 
understanding the causes of the differences between observed and 
measured response contributes to an understanding of the response of 
the ocean in the near surface layer. A third reason for studying such 
storms is that the deep ocean response may be typical of numerous 
storms of near hurricane force. If so, the associated changes in mass 
structure are significant at great depths. For example, the 
perturbation of the isopycnals will have drastic effects on acoustic 
propagation at depths at least as great as the SOFAR channel. The 
analysis also serves as an addition to a catalogue of techniques useful 
in analyzing transient oceanic behavior. 

Circular symmetry is assumed in the model used in this paper. 
However, it can be seen from Figure 1 that insufficient stations have 
been taken to positively identify the center of the vortex. Station 24 
was made after the vortex had been crossed once, when the research 
vessel returned to take a center station. It was felt that the station 
was taken close to the center but the possibility that it was not 
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cannot be excluded. The stationarity of Kara increases the amplitude 
of the response. Since the data were collected several days after the 
storm passed, it may be assumed that most of the energy is found in 
geos trophic motion. 



II > METHODOLOGY 



The existence of a cold water dome after the passage of a hurricane 
is suggestive of upwelling. This was concluded by Leipper (1967) and 
is supported by data collected by Landis (1966). Leipper^s data 
inspired O’Brien and Reid (1967) to pursue their study of upwelling in 
a two layer ocean as a response to a stationary hurricane. O’Brien 
and Reid established through a transient, non-linear model that the 
wind stress associated with storms of hurricane force and structure is 
capable of producing velocity divergence in the surface layer and large 
scale upwelling. O’Brien (1967) added turbulent mixing to the model 
and produced results which bear resemblance to Leipper ’s observations 
of the near surface layer. 

The problem can be approached somewhat differently. It is well 
known that surface divergence and associated upwelling will stretch 
water columns. A fluid will respond by spinning, thereby generating 
relative vorticity about the vertical axis. The effect extends to 
the bottom of the fluid and is transmitted much more efficiently than 
momentum transfer by friction. The model in this study utilizes the 
conservation of potential vorticity to describe the interior response. 
The assumption of a frictionless interior makes the continuously 
stratified problem tractable by allowing use of the conservation 
equation in the interior. This model and O’Brien’s differ because 
they investigate different aspects of the response to axisymmetric 
forcing. O’Brien solves the governing equations numerically as an 
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initial value problem in time because he is interested in the evolution 
of the response. By its very nature, his two layer model is incapable 
of obtaining the structure of the interior. However, it is successful 
in duplicating many aspects of the time dependent response. This 
study seeks to model the steady-state response of the ocean independently 
of the generating mechanism using a realistic density profile. 

The model used in this study was derived by Rooth (1970) and is 
presented fully in Appendix B. It makes use of the non-linear, time 
dependent equations of motion in deriving the vorticity equation. 
Hydrostatic equilibrium is assumed for the formation of the divergence 
equation. It is assumed that the data were collected a sufficient 
length of time after the storm passed to ensure that motion was 
essentially geos trophic. Combining the divergence equation and the 
vorticity equation results in a governing equation expressing the 
conservation of potential vorticity in a geostrophic velocity field. 



+ 



N^(Zo) 



= 

^2 



0 . 



II-l 



§ is the lifting of a density surface from its undisturbed state (called 
the matching level), f is the Coriolis parameter, and N^(Zq) is the 
square of the Brunt-Vaisala frequency of the undisturbed state. Zq is 
the geometric height in the undisturbed state. The observed mass 
distribution perturbed from the undisturbed state. The equation is 
the result of a transformation from cartesian coordinates with geometric 
height as the independent variable in the vertical to a "natural** 
coordinate system where density is the vertical coordinate. Figure 2 
illustrates the relationship between § and Zq. Geometric depth is the 
sum of Zq and 5* Formally, Z = Zq(p) + 5 (ZQ,r), where r is the radius. 
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Figure 2. The relation between density ( P) , undisturbed depth (Zq) and 
stretching (C) . 
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Not only is the numerical solution considerably simplified in this 
formulation, but solving in terms of the stretching is intuitively 
appealing. It is the stretching which generates the relative vorticity. 

The horizontal Laplacian of § transforms directly from cartesian 
to cylindrical coordinates. Equation II-l is therefore the governing 
equation in the latter system also. Let 

UZq,v) = R(r)Z(Zo) 

in order to separate variables, and II-l becomes 

Z”R + ZV^R = 0. II-2 



2 2 
For notational convenience, N (Zq) will be written as N^. 

Rearranging II-2 yields 

f2 . z" = -v2r 
T “ ' 

Since a function of Zq can be equal to a function of r only if they are 
each constants, then 



f2 ^ = -v2r = 1,2 

Wl' Z -TT ^ 



II -4 



Equation II -4 implies 

£i Z" - k^z . 0. II-5 

Equation II-5 must be solved numerically because of the inclusion of 
the Brunt-Vaisala frequency as a function of Zq, Note that the 
selection of sign for the separation constant implies that solutions 
are sought which are exponential in form in the vertical and harmonic 
in the horizontal. Equation II -4 implies 

V^R + k^R = o, II -6 



Writing the Laplacian in full, one obtains 
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R" + R'+ k^R = 0, II-7 

r 

assuming axisymmetr ic motion. Equation II-7 is a form of Bessel's 
equation and has the solution 

R = B Jo(kr) + C Nq (k^) . II -8 

Nq is a Neumann function or Bessel function of the second kind. Since 
the solution must be finite at r = 0, the constant C equals zero. 

With no loss of generality, the radial solution is normalized to unit 
amplitude at r = 0. The separation, then, is the following: 

5 §Q • II-9 

There is an infinity of solutions dependent on the value of k^, the 
arbitrary constant of separation called the wavenumber. In the appli- 
cation of equation II-9 to this study, the solution can be reduced to 
the response to a single, dominant wavenumber. The technique will be 
discussed following a summary of the numerical method used to solve 
equation II-l. 

Equation II-l was solved numerically on the IBM 370/155 computer 
at the University of Miami. A pre-programmed scientific subroutine 
package was utilized which is a modification by Hamming of Milne's 
predictor-corrector method. Since the predictor-corrector method is 
not self-s tarting , a fourth order Runge-Kutta starting procedure is 
used to generate sufficient points (Ralston, 1965). The numerical 
method requires specification of 5 and as initial conditions. 

At the bottom the stretching of an isopycnal is zero. It is assumed 
that a dominant wavenumber for the disturbance exists in the interior. 
It will be shown later that this assumption is realistic. There are, 
then, two arbitrary parameters, wavenumber and slope at the bottom, 
which have to be determined. A solution for specific values of slope 
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and wavenumber can be obtained by devising a scheme which matches the 
numerical solution to the density structure in two dimensions. The 
matching uses up the two degrees of freedom associated with the 
arbitrary parameters. A numerical scheme could be devised which would 
optimize the matching process by some best fit method. It was felt, 
however, that a subjective graphical method was justified considering 
the quality of the data. The approach chosen was to select an inter- 
mediate density surface and vary the input to the computer program 
until the amplitude of the response and the decay of the jQ(kr) Bessel 
function matched the observed perturbation. Essentially, the numerical 
solution, generated as an initial value problem, was forced to meet 
conditions on a surface on which density is conserved. 

The governing equation (II-l) requires the Brunt-Vaisala frequency 
of an undisturbed ocean as input. The input was taken from the 
temperature and salinity tables from R/V ATLANTIS stations 5446 and 
5447 (Fuglister, 1960). The stations were selected because they display 
temperature-salinity characteristics fairly typical of the Sargasso 
Sea and are relatively unperturbed by local variations. Figure 3 shows 
the observed response of the ocean to the storm. 
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Figure 3. Observed mass structure of ocean after passage of Hurricane 
Kara. Isopleths are values of sigma-t (<^ ^) • 
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Ill, DISCUSSION OF RESULTS 



A. DISCUSSION OF SOLUTION AND 
CONSISTENCY OF RESULTS 

Figure 4 shows a comparison between the computed and observed 
curves. The observed mass distribution has been shifted to a common 
vertical axis. This step is consistent with the assumption of 
axisymmetric motion implied in the transformation of the Laplacian 
in equation II-6 to cylindrical coordinates. As seen in Figure 3, the 
axis of the observed curves slants slightly away from the vertical. 

With the exception of the near surface layer and a few isolated spots 
at the sides, where the assumption of axisymmetric motion is in doubt, 
the difference between the curves is no more than 10% of the total 
perturbation amplitude (trough to crest) of any isopycnal. In most 
cases the percentage is considerably lower. The cross-hatched area 
in Figure 4 indicates the areas where the deviation of the two curves 
is significant. Considerable flattening of the density surfaces occurs 
at least as deep as 230 meters , and there is no evidence of upwelling 
above 100 meters . 

The lack of information about the state of the near surface layer 
prior to and during the storm precludes a complete analysis of the 
causes of the deviation between the two curves . According to Fuglister 
(1947) , average surface temperatures for October in the area of the 
storm range from 25.8°C to 26.1°C. These temperatures correspond very 
well to the surface temperatures found after passage of Kara. If the 
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CALCULATED 

OBSERVED 



Figure 4. Comparison between the observed mass structure and the computed 
response after Hurricane Kara. jQ(kr) = 0 at 80 km. 
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near surface layer did not cool significantly, then one must ask if 
upwelled water penetrated the thermocline and was mixed in the upper 
layer. It is possible that the apparent lack of surface cooling in 
the presence of obvious upwelling is the result of either horizontal 
advection in the surface layer moving the cold water away from the place 
of formation, or advection in the interior. The former is more likely, 
but the latter can partially explain the slanting of the axis of 
rotation of the vortex and the fact that the dome is located significantly 
to the south of the published hurricane track. 

Upwelling clearly took place as a response to Kara. One would 
expect to find evidence of downwelling at the sides to satisfy continuity 
requirements. O’Brien (1967) finds downwelling in a numerical model 
and Leipper (1967) finds evidence of it from the changes in shape and 
depth of the thermocline. Surface data taken after Kara do not indicate 
downwelling. The mixed layer depth is remarkably constant across the 
vortex and is typical of the ocean in this region. At the sides of the 
vortex, in the interior, downwelling definitely occurred in balance of 
the upwelling in the center region. This can be seen by the dipping 
of the isopycnals below the undisturbed levels Ths disparity between 
the curves at the top of the pycnocline suggests significant downward 
penetration of momentum below the depth of the mixed layer, leading to 
an effective Ekman depth of about 230 meters. That is, if the Ekman 
layer is considered that region where friction effects dominate, and 
where the conservation equations are not applicable, then the depth 
of significant differences between the observed and computed curves 
may be identified as the Ekman depth. 
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It is possible to check the consistency of the results with 

existing ideas of surface wind stress by working backwards from the 

theoretical displacement of the isopycnals. The technique begins with 

the Ekman pumping relation 

pfw = curl T , III-l 

where p is density, w is vertical velocity, and T is wind stress. 

Stokes* theorem may be written as 

Jt, dS = JJcurl T-dA, III-2 

S A 

or, by substituting equation III-l for curl t, as 

Jx-dS = pJJwf-dA. Ill -3 

S A 

Assuming an axisymmetric wind stress distribution, equation III-3 
becomes 

pwfnr^ = 2TTrT , III-4 



or 



T = wfr/2. 



Ill -5 



If w = 6§/6t where 6^ is the displacement of the most perturbed 
density surface and 6t is the time the storm was stationary, then the 



average stress distribution is 

T = 

2-6t 

The average stress may be found by taking an average over the curve 



R 2tt 

f = i • i r r gj-r -d 0 -dr 

2tt R II 2-6t 

= 1 • l2 f f ^ -r-d0-dr 

2tt R ^ ^ 2- 6t 



III-7 

III-8 



To simplify the integration, it was decided to utilize the series form 
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of the Bessel function: 
Jq(Z) = 1 



22(1!)Z 2^(2!)Z 26(31)2 



III-9 



Noting that 



and us ing 



1 

2tt 



2n 

! 



0 



1 , 



X = kr , dr = dX/k 

to change variables and integration limits, then III-8 becomes 

T = JQ(x)dX 

k"-R 0 2'8t 



III-IO 



The dominant wavenumber found by matching the theoretical and observed 
curves is 3.0 x 10"^. Therefore, the zero crossing occurs at r = 80 km. 
Using the well known relation 

^ = Pa'^* = PwV* = CdU2q . P3 III-ll 

enables one to use the average stress calculated in equation III-IO to 
find wind velocity, or 



u2 = 

10 Cp-Pa 



III-12 



where u is the friction velocity in air, v is the friction velocity 
in water, is the density of air (1.3 x 10*^ gra cm"^) , is the 
density of water (1 gm cm"^) , is the drag coefficient, and is 

the wind velocity at a specified height. The wind velocity is highly 
dependent on the choice of drag coefficient. Unfortunately, there are 
numerous estimates of the dependence of drag coefficients on wind speed 
and few agree. From Roll (1965) there seems to be a bunching of 
estimates for near hurricane force winds at Cq = 2.5 x 10"^. This is a 
figure 0*Brien (1967) also uses and terms "classical.” The results of 
this calculation for various values of the drag coefficient are 




• # 













tabulated below. Cj) = 1.5 x 10"^ is a widely quoted value which Malkus 
and Riehl (1960) found worked well in their model. 2.0 x 10"^ is the 
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mean of the two values. 

Cd Uio(KTS) 



1.5 


X 


10-3 


57.6 


2.0 


X 


10-3 


50.2 


2.5 


X 


10"3 


45.2 



Air Force penetrations of the storm during the time it was quasi- 
stationary found wind speeds between 45 and 55 knots. According to 
Colon (1966) , the wind distribution in the lower few thousand feet is 
uniform, or nearly so. Therefore the wind speeds should be repre- 
sentative of those near the surface. The tabulated results of the 
calculations are inconclusive. One may say there is no inconsistency 
between the theoretical result and gathered data. Occasionally the 
results of theoretical models are used to fix parameters upon which 
the models critically depend. A definite statement cannot be made in 
this case. However, the results indicate that the drag coefficient 
probably falls within the above range for winds of near hurricane force. 

As will be shown later in the paper, the penetration of the 
response of high wavenumber motion is less than that of lower wave- 
numbers due to the effect of the wavenumber on the exponential-like 
solution in the vertical (equation II-5) . The ocean acts as a low- 
pass filter to the response. This means that the difference between 
the two curves at the level where they diverge significantly may be 
due partially to the effect of the response at higher wavenumbers, 
which will change the height of a density surface in some unknown fashion. 
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It is also possible that the presumed center of the vortex was 
not the exact center. If not, then missing the center can account for 
some of the difference between the curves. 

B. THE MIXING PROCESS 
AND MOMENTUM TRANSFER 

Several calculations can be made which place an upper limit on 
the effects of the mixing process. The calculations assume a balance 
in the vertical between the downward diffusion of heat due to mechanical 
mixing and the upwelling of cold water from below the local thermo- 
cline. The governing equation for this process is 

w3t/ 3Z =K-32x/az2^ III-13 

where w is vertical velocity, T is temperature and < is the eddy 
diffusivity of heat. The solution of III-13 is 

T = To exp(w/K)Z. III-14 

</w is the e-folding distance for the temperature decay and is therefore 
the scale height of the local thermocline in this problem. H-j, the 
thermocline scale height, may be found by taking the logarithm of III-14, 
obtaining the following: 

2/rn: = 2jiTq + Z/Hx- III-15 

Plotting vs Z and finding the slope determines the scale height. 
Equation III-15 is the equation of a straight line on semi-log paper, 
and 1 /Hx is the slope of the line. Figure 5 is such a plot. Since 
pure mixing is assumed, the vertical velocity may be approximated by 
6h/6t, the difference in lifting between the observed and calculated 
curves divided by the relevant time scale of one day. < determined in 
this manner is 1400 cm^/sec. This value of k is several orders of 
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Figure 5. Log-linear plot of temperature (°C) vs depth (m) . 
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magnitude too large for abyssal processes but is probably the correct 
order for surface processes where eddy viscosity and eddy diffusivity 
may be nearly the same. If the point of divergence of the curves is 
at 230 meters, the eddy viscosity found from the Ekman depth formula, 
dQ = nV 2v/ f , is 2000 cm^/sec. 1400 cm^/sec is almost certainly an 
overestimate of the eddy diffusivity since Ekman divergence effects 
are neglected. If the depth of the mixed layer is neglected in the 
calculation of the eddy viscosity, dQ = 180 meters and v = 1200 cm^/sec. 
This calculation assumes the existence of two separate regimes: the 

mixed layer where eddy viscosity coefficients are very high; and the 
layer of Ekman divergence below the frictional layer, indicated by the 
disparity between observed and computed results. The similarity between 
the values of eddy viscosity and eddy diffusivity is presented for 
completeness and to suggest topics for further investigation. 

Under conditions of negligible advection, it is possible to 
extend an energy calculation by Turner (1969). Assuming pure mixing, 
the change in potential energy between the observed and theoretical 
curves can be used to work backwards to find wind velocity. Selecting 
an isopycnal where the two curves do not diverge and calling the density 
pQ a reference density, one can plot P^'PO vs (depth squared). 

Taking one half the area under the resulting curve solves the energy 
equation , 

EpoT “ gJCPZ’Po^^ ^ III -16 

graphically. Figure 6 is a plot of Ap vs Z^ for the center station, 
where maximum upwelling occurred. The change in potential energy found 
planimetrically is 10.14 x 10^ ergs cm"^. By dimensional arguments. 
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Figure 6. Density difference (gm/cm^) plotted vs square of depth (cm^) 
for potential energy calculation. Station is compared to the center 

of the computed curves . 
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Turner asserts that it is reasonable to write the time rate of change 
of energy as 

= C^TUio = C2pUio^- III-17 

The tabulation below gives values of wind speed for two values of the 
constant C 2 found by Turner. 

C2 Uio(KTS) 



2.6 X 


10"8 


71.4 


1.3 X 


1-* 

0 

1 

00 


90 



Both values of wind speed are too high if pure mixing is responsible 
for the potential energy change. If one accepts 50 knots as a 
reasonable value of wind speed during the storm, then surface mixing 
accounts for about half the energy change. That is, at least half of 
the density field discrepancy from the theoretical curves must derive 
from the deep Ekman divergence. As the drag coefficient. Turner’s 
constants are not well defined. The dependence of Turner’s constants 
and the drag coefficient oxi wind speed is unknown. Turner’s method 
requires a one -dimensional process and a regime in which surface cooling 
and resultant convective overturning are unimportant. Both processes 
are usually significant under hurricane conditions, although their 
effect may be minimal in this case since surface cooling apparently did 
not occur. Selecting the center station minimizes the advective effects, 
since § is a maximum there and 55/^r approaches zero. As Turner 
correctly points out, his method is illustrative but hardly conclusive. 

C. THE EFFECT OF LENGTH SCALE ON RESPONSE 



Some interesting properties of the solution as a function of scale 
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of forcing and depth are illustrated in Figure 7. The governing equation 
(II-l) was solved using the same initial conditions and density profile 
but different values of wavenumber. The length scale was determined 
by selecting the value of radius which corresponded to the first zero 
of the Jq Bessel function, found from the formula (1 /wavenumber) x 2.4 « 
The values of 5 were normalized to 100 meters of Ekman pumping at 
Zq = 25 meters and isopleths of constant § were plotted. The response 
may be split into a barotropic and baroclinic regime. In the barotropic 
regime, response is characterized by the stretching being a linear 
function of depth, since stretching is proportional to vertical velocity. 

The transition from a baroclinic to a barotropic response is rather 
broad in scale, occurring between 300 to 500 km. For scales greater 
than 500 km, the response may be continued to be barotropic. The 
response for high wave numbers is confined to a shallow surface layer. 

The latter provides the rationale for basing the solution of the 
vorticity equation on the existence of a dominant wavenumber in the 
interior. The dependence of the vertical structure of the response on 
the scale of the motion suggests a definite problem in the representation 
of a continuously stratified ocean by a two-- layer model. For the higher 
wavenumbers, the two layer model lacks sufficient resolution to describe 
the response adequately. 

The ordinate in Figure 7 is undisturbed depth (Zq) . Since depth 
is equal to the sum of the undisturbed depth and the stretching, zero 
depth in the figure should not be confused with the free surface. The 
free surface lies considerably below Zq = 0. 

It is anticipated that similar diagrams can be generated using 
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Figure 7. Isopleths of stretching, 5 (m) , as a function of depth and 
length scale. Values of ^ are normalized to 100 meters pumping at 
25 meters matching depth. 
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density distributions typical of other areas of the ocean. The diagrams 
should prove useful in predicting the response in these areas. A 
density profile with higher resolution than that obtained from 
Fuglister's Atlantic Ocean Atlas is preferable. 



D. GEOSTROPHIC VELOCITIES 



It is possible to calculate geostrophic velocities using the 
solution of equation II-l. Noting 

Z(p,r) = ZqCp) +§(Zo,r), III-20 

and from Appendix A , 



5^ _~5z\ 5p 

3rjzo 3rjp ’ 

one may use the thermal wind equation, 

M 

hZQ pf hrjzQ 

to find velocities throughout the column and for any radius. 



Substituting III-21 


into 


ni-22, . 


one 


obtains 






= _&|z\ 








^0 


pf ^jp 




^0 


Equation III-23 may 


be integrated 


in 


the vertical: 




’Jzo 


Zq 

du = i r 




!e • f dz 
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f J,. 
^0 


P c 


)zo yirjp 



III-21 



III-22 



III-23 



III -24 



Noting (bottom) = 0 and N^ = , equation III -24 becomes 

0 G p qz Q 



Uz(, = -i nJ dZo 



III-25 



Introduce III-20 for Z in III-25, obtaining 



Uz' = -- Nq dZQ . 

^0 f i" 0 \hr)p ° 



III-26 



Geostrophic velocities were calculated numerically using equation III-26 
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and the result established previously, 5 = |q • The derivative 

of I is 



3§/ar = k §0 Ji(kr) . III-27 

Equation III-26 is a valid expression for the velocities induced by 
Kara, except in the Ekman layer, where the predicted density field 
departs from the observed. Veronis (1956) showed that motion which 
has forcing on scales larger than the baroclinic radius of deformation 
is essentially geostrophic. In this case, the scale of forcing, 80 km, 
is more than twice the baroclinic radius of deformation, 36 km. 

Figure 8 is a comparison of the geostrophic velocities based on 
equation III-26 and velocities calculated using the dynamic method 
(called the observed velocities) . The reference level for calculating 
the observed velocities was found using Defant's method. The level 
was conspicuously present in all the profiles. Although the percentage 
of error is high at the lower depths, the absolute magnitude difference 
is low. In general, the magnitude of error remains the same as depth 
decreases. The two velocity profiles agree rather well in magnitude 
and general shape. Aside from the reasons for the difference between 
the curves in the Ekman layer, which have been discussed, the calculation 
of geostrophic velocities poses special problems. In this case, the 
distances between stations used to calculate the velocities by the 
dynamic method were extreme. It was found that velocities at a given 
radius and depth varied as much as 5 cm/sec depending on the particular 
combination of stations used. The curves represent the smoothest 
combination. The observed velocities were computed from the non-centered 
observed data, which certainly accounts for some of the difference 
between the curves. The effect of centripetal acceleration has been 
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Figure 8. Comparison of numerically calculated geostrophic velocities 
(cm/sec) with velocities found using the dynamic method (observed). 
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neglected in the calculation of velocities, and in the formulation of 
the vorticity conservation model. The inclusion of centripetal 
accelerations makes the model non-linear, introducing an undesired 
degree of difficulty in solving the vorticity equation. A measure of 
the importance of the centrifugal force term is the Rossby number, 

f *r 

O^eumann and Pierson, 1966). Choosing r = 40 km, v = 40 cm/sec, and 
f = 7.5 X 10"5 sec"^, the Rossby number is 0.13. It can be shown that 
geos trophic velocity is an overestimate of the gradient velocity if 
curvature is important (Haltiner and Martin, 1957). 

The success in obtaining substantial agreement between the 
theoretical mass field and velocities and the observations supports the 
validity of the initial model formulation. 



IV , SUMMARY AND CONCLUSIONS 



A model based on the conservation of potential vorticity was 
solved numerically. The solution was utilized to investigate some 
aspects of the oceanic response to a tropical storm as a function of 
horizontal scale. In particular, the interior response of a moderately 
intense, quasi-stationary cyclone was duplicated. The matching was 
successful below a depth of 230 meters which was associated with the 
Ekman depth. The lack of information concerning the state of the 
ocean in the near-surface layer prior to the passage of Hurricane Kara 
compounded the problem of discovering what caused the divergence of 
the observed and computed curves. The divergence was attributed to the 
neglect of high wavenumber modes when solving for the interior response , 
and to some ill-defined combination of vertical mixing and momentum 
transfer by friction. Gross estimates indicate that the two processes 
are of equal importance in this particular case. An upper limit for 
the eddy diffusivity is 1400 cm /sec. Vigorous mixing, if it occurred, 
may not have resulted in significant cooling of the mixed layer. 

Geostrophic velocities were calculated from the raw data using 
Defant*s method to find a reference level of no motion. These 
velocities were compared to those calculated from the numerical solution. 
The velocities agree rather well below the near-surface layer. Since 
the velocity computation integrates the stretching, the comparison of 
velocities should not be expected to be as successful as the comparison 
of I itself. The integration of § means that errors are accumulated 
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over the entire column* Geostrophic vertical shear matches better 
than geostrophic velocities. It is plausible that the nearly constant 
error in magnitude between the velocity fields is a result of 
inaccuracies in the dynamic computations; specifically, limitations in 
the Defant method of obtaining the level of no motion may be the cause. 

It was found that the scale of forcing can have dramatic effects 
on the penetration of the response. For equal Ekman pumping, higher 
wavenumber (short scale) response does not penetrate very deep into 
the ocean compared to low wavenumber response. Thus the scale of the 
disturbance must be considered when choosing the depth and number of 
layers for modelling a continuously stratified fluid. A two layer 
model may lack the resolution necessary to accurately describe the 
response in the upper layer. Adequate modelling of the near-surface 
layer by a model such as the one used in this study requires the 
inclusion of high wavenumber motion. 



APPENDIX A 



THE NATURAL COORDINATE SYSTEM 
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Although the standard cartesian coordinate system is the one 
normally used in hydrodynamics , there is no reason why the vertical 
coordinate must be geometric height, a fact recognized long ago by 
meteorologists and some oceanographers. The use of a "natural” vertical 
coordinate, such as pressure, density or potential density, can simplify 
an analysis but render boundary conditions formidable. 

In the potential vorticity conservation model, density will be 
utilized as the "natural" vertical coordinate. This is advantageous 
because it enables one to solve directly for the displacement of an 
isopycnal from equilibrium. The method of coordinate transformation is 
discussed in Haltiner and Martin (1957). It is summarized below since 
the method is implied in the model and is used directly in deriving the 
geostrophic velocities (equation III-26) . 

Consider 0(x,y,z,t) = constant. The change of Tj(x,y,z,t) along 
0 = constant is equal to the change along the horizontal plus the 
change along the vertical due to the slope of 0, i.e. From Figure 9, 




A-1 



Now , 




A-2 



by the chain rule at a constant time, tQ. A-2 may be written as 




A-3 



where is the slope of the 0 surface in the x direction. Let T] - p 



3x 

and A-3 becomes 










I 



A-4 
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Figure 9. 0 = constant surface separated into horizontal and vertical 

components to illustrate coordinate transformation technique. 
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APPENDIX B 

THE POTENTIAL VORTICITY CONSERVATION MODEL 
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The governing equation for the interior motion of a hurricane, it 
is hypothesized, expresses the conservation of potential vorticity. 

The equations of motion used are non-linear, time dependent, hydrostatic 
equations in a rotating system characterized by the Coriolis parameter, 
f, which is a function of the horizontal coordinates alone: 





u • 


Vu 


- fv 
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0 

II 


B-1 
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U • 


W 


+ fu 


+ 


II 

0 


B-2 



where tt = ^ + gZ . The vorticity equation may be formed by cross- 
differentiating equations B-1 and B-2, resulting in B-3: 

fe ^ I 7 ] + vfy = 0 B-3 

B-3 may be written as 

[It * “ I? y 

since f is not a function of time. That is, 

^ (f4v^-Uy)+(f4vx-Uy)(u3j4Vy) = 0 , B-5 

where 



d _ a . a . a 

— r: + u + V ^ . 

dt “ ot Ok ay 

' a 

The divergence equation may be formed by taking (equation B-4) 



and adding to ^ (equation B-2) . 

ay 



^ — (u , +7 ) -f (v -u ) + V^TT + f ^u-f v+u^“+v^+v u +u V “0 B -6 
y^ y X X y X y y x 

Assuming geostrophic motion reduces the balance equation to 



f(v^-Uy) = V^TT. B-7 

The continuity equation is 

u +v +w =0. B-8 

2^. y ^ 

In the vertical, w depends only on z(p), therefore 
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+ ^ • ^ = 0. 

^ 3p 3z 

Noting w = dz/dt, equation B-8 becomes 
^(u^4v ) + = 0. 



B-9 



Substitute equation B-9 into the vorticity equation (B-4) to obtain 
, d 



^(Sz/3p) 

- u ) • 3t_ = 0. 

y 3z/5p 

Rearranging terms, equation B-10 becomes 

^(f + Vx - Uy) _ ^(3z/ap) ^ ^ 



f + v„ - 






3z/3p 



or, by definition of the logarithm, 

^Vy-U^T+f 



dt\az/3p 



B-ll 



B-12 



Equation B-12 implies that for each fluid element 

= constant. B-13 

3z/5p 

The existence of a uniform reference state implies that the constant is 

a function of density alone. If the reference state is rest, relative 

vorticity is zero, and the constant, 0(p) , equals f/(3z/5p)Q. The 

vorticity equation, then, becomes 

^x-^y+f = _f . B-14 

3z/3p (9 z/3p)q 

Substituting equation B-7 into B-12 and rearranging, one obtains 

fi+^^TT VazN = ^ ^ b-15 

\ f2 Aap/0 ap 

Next, apply the hydrostatic equation, 

TT = + gz • B-16 

p ^ 

Rearranging and taking the partial derivative with respect to z results 



in 



|^(TT-gz)p = P^. 



B-17 
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Assuming hydrostatic equilibrium, B-17 becomes 



dz 



.(TT-gz)p = -gp. 



B-18 



Integrating B-18 results in the following: 

Z 

prr-pgz = PQ ^0 “ J PSdz. 

0 

The last term in equation B-19 can be integrated by parts: 

Z p 

-Jpgdz = -pgz + Jgzdp . 

0 PO 

Substituting B-20 into B-19, the desired forms of the hydrostatic 

equation is achieved: 

^ = ^0 + T fs2dp . 

^ PO 

Take the Laplacian of B-21: 

p 

V^TT = + 1 JgV^zdp . 

Po 

Substitute equation B-22 into equation B-15: 

P 



B-19 



B-20 



B-21 



B-22 



h. 1 (v2tt + i Jgv2zdp)(^|£\ 

l,3pjo ° '>io V5'>7o 



PO 

Therefore, if Pq is chosen as the surface density, where h is the free 
surface height, 



h. .(i£\ . S^(v2h + i fv2zdp)(|^)„ . 
ap UJo ppo 



B-23 



The depth, z, may be split into two parts, 

Z = Z(j[p) + |(x,y ,Zq) . 

Equation B-23 becomes, then, after rearranging, 

v2h + - fv2|dp = |i . 

pi. g 9o \^/0 



Po 



B-24 



Next, take the derivative of B-24 with respect to p: 



- Ofl (^' 
p gpo yap 



(1^0 ■ 
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B-25 
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5p 
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SO equation B-25 becomes 

Po PoN2(zo) azo2 ’ 

or 



7^5 + 



f2 

N^(Zq) 
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Equation B-26 is the governing equation used in the analysis described 
in the text. The basic derivation is attributed to Rooth (1970). 
Several steps were expanded by this author for clarity and completeness 
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